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Abstract. We revisit the kink-like parametrization of the deceleration parameter q(z) [ 1 ], 
which considers a transition, at redshift zt, from cosmic deceleration to acceleration. In this 
parametrization the initial, at z^> Zt-, value of the q-parameter is q^ its final, z = —1, value is 
qj and the duration of the transition is parametrized by r. By assuming a flat space geometry 
we obtain constraints on the free parameters of the model using recent data from type la su¬ 
pernovae (SN la), baryon acoustic oscillations (BAO), cosmic microwave background (CMB) 
and the Hubble parameter H(z). The use of H(z) data introduces an explicit dependence of 
the combined likelihood on the present value of the Hubble parameter iTo, allowing us to ex¬ 
plore the influence of different priors when marginalizing over this parameter. We also study 
the importance of the CMB information in the results by considering data from WMAP7, 
WMAP9 (Wilkinson Microwave Anisotropy Probe - 7 and 9 years) and Planck 2015. We 
show that the contours and best fit do not depend much on the different CMB data used and 
that the considered new BAO data is responsible for most of the improvement in the results. 
Assuming a flat space geometry, qi — 1/2 and expressing the present value of the decelera¬ 
tion parameter qo as a function of the other three free parameters, we obtain zt = 0.67^QQg, 
r = 0.26 ^q'^q and qo — ~0.48^q'^ 3 , at 68 % of confidence level, with an uniform prior over Hq. 
If in addition we fix qj — — 1 , as in flat ACDM, DGP and Chaplygin quartessence that are 
special models described by our parametrization, we get zt — 0 . 66 ^q'q 4 , t = 0 . 33 ^q'q 4 and 
go = -0.54l°;°7> in excellent agreement with flat ACDM for which r — 1/3. We also obtain 
for flat reCDM, another dark energy model described by our parametrization, the constraint 
on the equation of state parameter —1.22 < w < —0.78 at more than 99% confidence level. 

Keywords: cosmological parameters - cosmology: supernova type la, baryon acoustic oscil¬ 
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1 Introduction 

Since the end of last century, when type la supernovae (SN la) observations [2] have suggested 
that the expansion of the universe is speeding up, we have seen a growing number of evidences 
confirming this cosmic acceleration. In fact, much progress has been done since then, but 
understanding the fundamental physics behind cosmic acceleration is still an open question 
and one of the main challenges of modern physics. 

Several experiments are underway and being planned and many theoretical approaches 
have been suggested to investigate cosmic acceleration. In this work we follow the so called 
“kinematical approach” [1, 3, 4] which is based on a parametrization of the deceleration 
parameter q as a function of the redshift z. The advantage of using a g(z)-parametrization 
instead of, for instance, a parametrization for the equation of state, is due to its generality. 
In this case, we are not restricted to general relativity and the assumptions regarding dark 
matter and energy are minimum. We only have to assume a metric theory of gravity. Natu¬ 
rally, we pursue a parametrization which can describe several cosmological models in a wide 
redshift range. 

Here, our main interest is in the transition from cosmic deceleration to acceleration. 
We are particularly interested in the redshift of this transition, its duration (in redshift) 
and the present value of the deceleration parameter. To achieve our goal we chose the 
kink-like q(z) expression proposed in [1], which is quite general. For the sake of simplicity 
we assume a flat FLRW metric. To constrain the parameters of the model we use the 
following kinematical tests: SN la, BAO/CMB and also observational Hubble data (OHD). 
The BAO/CMB test is based on the ratio between the comoving angular-diameter distance 
to the photon-decoupling redshift g?a(^*) and the “dilaton scale” Dv(zbao)- It is equal to 
zb AO times the ratio 7 Z(z*) / A(zbao)-> where 1Z is the CMB shift parameter and A is the BAO 
scale [31]. The OHD test was not used before to constrain our kink-like (/(z)-parametrization. 
With the inclusion of H(z) data we also investigate the influence, on the constraints of the 
parameters, of different priors (two gaussian and one uniform) on the present value of the 
Hubble parameter Hq. 
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This paper is organized as follows: in Section 2, we present the kink-like g-parametrization 
and discuss some of its properties. In Section 3, the cosmological tests (SN la, BAO/CMB 
and H(z)) we use to constrain the model parameters are presented. Finally in Section 4 our 
main results are exhibited and discussed. 


2 The model 


We investigate in this work flat cosmological models which deceleration parameter (q = 
— p = ^ (-^) — 1), after radiation domination, has the following redshift dependence [1, 4] 


q(z) = q f + 


(gi ~ Qf) 

<k f 1 + zt 
Qf V 1 + * 


1/ T * 


( 2 . 1 ) 


In the above expression zt is the transition redshift (q(zt) = 0) from cosmic deceleration 
(q > 0) to acceleration (q < 0), q\ > 0 and qf < 0 are, respectively, the initial (z » Zt) 

and final (z —>■ —1) asymptotic values of the deceleration parameter. The parameter r is 

2 • • • 

associated with the width of the transition and is related to the jerk (j := g fjr) at the 
transition 


T 


-1 



1) 

dq(z ) 

g // 

_dln(l + z) _ 




( 2 . 2 ) 


As usual, a = 1/(1 + z) is the scale factor of the Friedman-Robertson-Walker metric, Hq = 
100 h km/(s Mpc) is the present value of the Hubble parameter (H = and a dot over a 
quantity denotes differentiation with respect to the cosmic time. 

From the definitions of q and H we can write 


TJ U r 1 + gO®) A 

H =- Hq exp / —--ax, 

Jo 1 + x 


(2.3) 


that can be integrated by using (2.1) to give 


AMI 2 

. H 0 _ 


(1 + z) 


2(1+qi) 




(1 + Zt 

V 1 + z 


1/t 


Qi (1 + z t) l H 


-g/ 

g/ 


2 r(qi-qf) 


(2.4) 


We can relate the final value of the deceleration parameter, qf := q(z = —1), to its 
current value qo := q(z = 0) 


g/ 


Qi (1 + z t) l H 



1 - (1 + Zt) 1/T 


5 


and by substituting the above equation in (2.1) we get 


qoqi 

1 1 

T—1 

1 

1_1 


(1 +z)~ 1 / T 

qo + (i + Z t )r • 

- Qi 

- qo 


(2.5) 


(2.6) 
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As mentioned above, with our kink-like q(z) parametrization we aim to describe the 
transition from a cosmic decelerated phase to an accelerated one. Therefore, since qi > 0 and 
qf < 0, by using Eq. (2.5) it is straightforward to show that the parameter r is constrained 
to the interval 


0 < r < 


ln(l + z t ) 


In 


i-® 

qi 


(2.7) 


In most cosmological scenarios, large scale structure formation requires that at early 
times the universe passes through a matter dominated era in which, H 2 oc (1 + z) 3 and, 
consequently, q = 1/2. Since the g-parametrization we are working with is designed to 
describe cosmic evolution starting from a matter dominated decelerated phase on we fix 
qi — 1/2 in this work. In fact, with this assumption we do not lose much generality, and are 
essentially excluding from our description models with a coupling in the dark sector (see [1] 
for more on this point). We should remark that although the validity of our parametrization 
can be extended to high redshift (z <1-2 x 10 3 ) it is not valid during the radiation 
dominated era (RDE) when q — 1. However when dealing, for instance, with BAO and CMB 
data, we are implicitly assuming, besides the existence of baryons, a regular RDE and that all 
standard assumptions for the very early universe are valid. Finally, we should also mention 
that several flat cosmological models investigated in the literature are special cases of the 
kink-like q(z) parametrization given in Eq. (2.1), as for instance: w CDM, ACDM, the DGP 
brane-world model [5], the quartessence Chaplygin model [6-9] and the Modified Polytropic 
Cardassian model [10, 11] (see references [1] and [4] for more on this point). 


3 Observational Data 

In this work we derive the constraints set by three different experiments upon the parameters 
Zt , r and go- We use the following datasets: a) observations of type la supernovae from the 
Joint Lightcurve Analysis (JLA) of the Sloan Digital Sky Survey (SDSS-II) and the Supernova 
Legacy Survey (SNLS) samples [12]; b) observations of baryon acoustic oscillations in the 
SDSS DR7 [13] and SDSS DR9 [14], in the WiggleZ Survey [15] and in the 6dF Galaxy Survey 
(6dFGRS) [16]; c) measurements of cosmic microwave background temperature anisotropy 
from 2015 Planck [17] and WMAP9-year [18]; d) observations of the Hubble rate H(z) 
from [19-24]. In the following subsections, the statistical analysis used to deal with those 
observables is briefly described. 

3.1 Type la supernovae 

The kink-like g-parametrization was used in [1] to investigate discrepancies, reported in the 
literature, between constraints from different SN la datasets in a more general context. In [4] 
we used only one SN la dataset, the 288 SN la compilation from Kessler et al. [25] (sample 
“e” in their paper), analyzed with two light-curve fitters, MLCS2k2 [26] and SALT2 [27], 
and showed that the tension appeared to be between the fitters instead of the datasets. It 
has been argued in the literature (see [28] and references therein) that the origins of the 
discrepancy between the fitters are the rest frame U band calibration and the priors in SN la 
colors, and the community seems to agree that SALT2 is more reliable in these aspects. In 
this work we use the 740 SN la compilation from the JLA [12], available solely for SALT2, 
and follow the analysis proposed in Appendix E of [12], which we are going to briefly describe 
in this section. 
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A key quantity in SN la investigation is the distance modulus 


Hth(z; M) = 5 log [. D l (z\ 6)\ + Ho(h). 


(3.1) 


Here h is the Hubble constant in units of 100 km/s Mpc , 6 denotes the set of cosmological 
parameters of interest other than h, 


Ho(h) = 5 log 


^ 10 3 c/ (km/s) ^ 


= 42.38 - 5 log h. 


(3.2) 


and 


D L (z-h,e) 



dz' 

H(z'-J h 6)/Ho ’ 


(3.3) 


is the dimensionless luminosity distance (in units of the Hubble distance today), assuming 
flat geometry. 

In the context of SALT2, the j-th SN la distance modulus is modeled by 


Hj = m* B - - M b + axij - /3cj , 


(3.4) 


where the parameter c is a measurement of the SN la color, x\ is related to the stretch 
of the light curve and m* B is the peak rest-frame magnitude in the B band. Its statistical 
uncertainty is given by 


/ / r 
«) = K- 


J ) 2 + a 2 « 1 ) 2 + / 3 2 K ) 2 + 2 a (. 


rripXi > 


a- 


- 2 


-2a/3(af c ) 


(3.5) 


Here the quantities (cr™ 73 ) 2 , (crj 1 ) 2 , (crp 2 , g™ bC and crJ lC are the components of the 

covariance matrix of (m^, aq, c). We remark that m^, x\ and c are derived from the fit to the 
light curves while the global parameters cq /? and the absolute magnitude M are estimated 
simultaneously with the cosmological parameters and are marginalized over when obtaining 
the cosmological parameter constraints. 

In order to take into account the dependence of the absolute magnitude with host galaxy 
properties, it was introduced an additional parameter AM (originally proposed by [29]) such 
that 


M b = M b if M ste u ar < 10 10 M o , 

Mb = M b + A m otherwise, (3-6) 

where M ste u ar is the host galaxy stellar mass. 

The likelihood for SN la data is given by 

LsN= vmmw eM - xiN/2) ' <3 - 7) 

where C is the N x N covariance matrix of the N SN la observations. 

In this work we assume that the estimate of distances is independent of the estimate of 
cosmological parameters [12]. Under this assumption, the distance modulus can be approxi¬ 
mately written as 

fi(z) = (1 - a)lib + ocfib+1, (3.8) 

where a = \og(z/zt,)/ log( 25+1 /zj } ) and ///, is the distance modulus at 25 . 
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This function was fitted to the measured Hubble diagram by minimizing the y 2 function 

Xsn = (M- /^)±“V- p(z)). (3.9) 

The parameters to be determined are cq /?, A m and the values of /i^ in each bin (they chose 
31 log-spaced points in the redshift range 0.01 < z < 1.3). In order to obtain unique values 
for fi\) it was used a fiducial value of Mg = —19.5. 

Once the correlation matrix C for the parameters is obtained we can use it in the 
cosmological analysis by inserting in Eq. (3.7) the following function 

Xsn = (fJ'b-M - 5 log - M - 51og D L (z b ,0)), (3.10) 

where A4 = M — 5 log h+ 42.38. The cosmological constraints are obtained marginalizing over 
the parameter M with an uniform prior. We remark that in our analysis the full covariance 
matrices, including systematic errors, as reported by [ 12 ], were used. 

3.2 Baryon acoustic oscillations and cosmic microwave background 

In this section we are going to briefly review the combined BAO/CMB test. However, before 
doing so, some remarks are required. We emphasise that our approach is kinematical in 
which we avoid strong assumptions about the dark sector. All the tests we consider here are 
essentially based only on distances and we are not considering tests directly based on the 
growth of perturbations. Furthermore, as mentioned before, although the parametrization 
given in Eq. (2.1) is not assumed to be valid during RDE, we assume, besides the existence 
of baryons, a regular RDE and that all standard assumptions for the very early universe are 
valid. Thus, we can define the comoving sound horizon at the photon-decoupling epoch by 

. . c /‘ 1 /(i+ 2 *) da . . 

r s (z*) = / - =, (3.11) 

Vs Jo a 2 tf(a)yi + (3ft 60 /4tt 7 o)a 

where fl 7 o and D 50 are, respectively, the present value of the photon and baryon density 
parameter and z* is the redshift of photon decoupling. To obtain the BAO/CMB constraints 
another relevant quantity is the redshift of the drag epoch ( 2 ^), when the photon pressure is 
no longer able to avoid gravitational instability of the baryons. The WMAP9 [18] values for 
these two redshifts are z* = 1090.97^0 86 an d z d — 1020.7 ± 1.1, while the 2015 results for 
Planck [17] are z* = 1090.00±0.29 and z^, = 1059.62±0.31. Also important are the “acoustic 
scale” l a — 7 r d r A ^*) ’ an dthe “dilation scale”, Dy{z) := [d\(z)cz/H(z )] 1 ^ 3 , introduced in [31]. 
Here, gU(^*) = c fj* dz'/H(z') is the comoving angular-diameter distance. 

We are using the results from Padmanabhan et al. [13], which reported a measure¬ 
ment of r s /Dy at z = 0.35, (r s (zd)/Dy( 0.35) = 0.1126 ± 0 . 0022 ), Anderson et al. [14], 
(r s (z d )/D v ( 0.57) = 0.0732 ± 0.0012), Beutler et al. [16], (r s (z d )/D v ( 0.106) = 0.336 ± 
0.015) and Blake et al. [15], which obtained results at 2: = 0.44, z = 0.6 and 2: = 0.73 
(r s (z d )/D v (0.U) = 0.0916±0.0071, r s (z d )/D v ( 0.6) = 0.0726±0.0034 and r s (z d )/D v (0.73) = 
0.0592 ±0.0032). 

We combine theses results with the Planck 2015 [17] values 1000* = 1.04106 ± 0.00031, 
which gives l a = tt/0* = 301.77±0.09, r s (zd) = 147.41 ±0.30 and r 5 (z*) = 144.71 ±0.31, for 
the combined analysis TT, TE, EE±lowP±lensing (see Table 4 of that reference). Following 
[32], we multiply the value of r s (zd ) by a constant factor of 1.0275 to take into account the 
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difference between CAMB [33] results and the Eisenstein and Hu approximation [34] used in 
the BAO analysis. With this information we can write the y 2 for the BAO/CMB as 


XbaO/CMB — Xt C 1 -X-> 


where 


and 


X = 






dAz*) 
£V(0.106) 
dAz*) 
D v { 0.35) 
Aa{z *) 
D v ( 0.57) 
dAz*) 
D v { 0.44) 
dA{z *) 
D v ( 0.6) 
dAz*) 
D v ( 0.73) 


— 30.84^ 

— 10.33 
-6.72 
-8.41 

— 6.66 

- 5.43 , 


C" 1 


/ 0.52552 
-0.03548 
-0.07733 
-0.00167 
-0.00532 
-0.00590 


-0.03548 

24.97066 

-1.25461 

-0.02704 

-0.08633 

-0.09579 


-0.07733 

-1.25461 

82.92948 

-0.05895 

-0.18819 

-0.20881 


-0.00167 

-0.02704 

-0.05895 

2.91150 

-2.98873 

1.43206 


-0.00532 

-0.08633 

-0.18819 

-2.98873 

15.96834 

-7.70636 


-0.00590 \ 
-0.09579 
-0.20881 
1.43206 
-7.70636 
15.28135 / 


(3.12) 


(3.13) 


(3.14) 


is the inverse covariance matrix derived by using the equations given above together with the 
correlation coefficients r = 0.39, r = 0.453 and r = —2.12417 x 10 -7 calculated for the r s /Dv 
pair of measurements at z = (0.44,0.6), z = (0.6,0.73) and z = (0.44,0.73) as presented in 
[35]. In Ref. [4] the correlation between the measurements at z = 0.44 and 0.73 was assumed 
to be zero since it was not available in the original reference [15]. 

We also used the results from WMAP9-year [18]: l a = 302.35±0.65, r s (zd) = 152.3±1.3 
and r s (z*) = 145.8 ± 1.2 which gives 


X = 






dA(z *) 

ZV(0.106) 

dAz*) 

D v { 0.35) 
dAz*) 
D v { 0.57) 
dAz*) 
D v (0.U) 

dAz*} 

.Dy (0.6) 
dA{z*) 
D v ( 0.73) 


- 30.96 ^ 

- 10.37 
-6.74 
-8.44 
-6.69 

- 5.45 , 


(3.15) 
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and 


C- 1 


0.50557 

-0.29215 

-0.63687 

-0.01373 

-0.04382 

-0.04862 


-0.29215 

20.60694 

-10.33210 

-0.22271 

-0.71095 

-0.78883 


-0.63687 

-10.33210 

62.47146 

-0.48551 

-1.54983 

-1.71962 


-0.01373 

-0.22271 

-0.48551 

2.87954 

-2.99476 

1.38827 


-0.04862 \ 
-0.78883 
-1.71962 
1.38827 
-7.75021 
15.04645 / 


-0.04382 

-0.71095 

-1.54983 

-2.99476 

15.74975 

-7.75021 


3.3 H(z) 


(3.16) 


Constraining cosmological models with observational Hubble data (OHD), i.e., by direct 
measurements of the expansion rate, has received considerable attention in recent years. The 
two main methods to measure the Hubble parameter as a function of redshift, H(z), are the 
differential age method and the radial BAO size method. For review, including a discussion 
on systematic uncertainties and possible source of errors, see [36]. 

The differential age method has been proposed in [37]. It has been shown in this 
reference, that by measuring the age differences between two ensembles of passively evolving 
galaxies at somewhat different redshifts, it is possible to obtain the derivative dz/dt and this 
quantity is directly related to H(z), 


H(z) = 


1 dz 
1 ~\~ z dt 


(3.17) 


By using the radial BAO peak position as a standard ruler in the radial direction, in Ref. 
[20], for the first time, it has been obtained a direct measurement of H(z) by the second 
method mentioned above. Recently, Farooq and Ratra [38] compiled a list of 28 independent 
measurements of the Hubble parameter at different redshifts and used these measurements 
to constrain cosmological models. Here we use their compilation (see Table 1 of [38]), but 
substitute the result from [39] (at z = 2.3) with that from Delubac et al. [23] (at z = 2.34). In 
[23] it was used over 137500 quasars from Data Release 11 of the SDSS-III Baryon Oscillation 
Spectroscopic Survey (BOSS), more than twice the number of quasars considered in [39]. 
This is the most important improvement of [23] over [39]. Other treatment improvements 
are described in [23]. To avoid possible correlations and to assure data independence, we also 
excluded the three data points from [40], at z = 0.44, z = 0.60 and z = 0.68, since the same 
data is used in our CMB/BAO analysis [15]. In our compilation we have also included two 
new measurements from cosmic chronometers at redshifts 1.363 and 1.965 from [24]. 

In order to constrain the 6 parameters of the model we use the following \ 2 h function 


27 

xh = T 
2=1 

where Hth(zi,Ho,0) denotes the Hubble parameter at Z{ predicted by the models, H 0 b s ^ is 
the i-th measured one and an,i is its uncertainty. 

In this work we marginalize \ 2 h over (see the Appendix for details), using three 
different priors: 

1. an uniform prior over all Hq possible values. 

2. A Gaussian prior using the measurement from Riess et al. [41], Hq = 73.8 ± 2.4 km 
s -1 Mpc -1 ; 


H th (zi,Ho,e)-H t 


obs,i 


H.i 


(3.18) 
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3. A Gaussian prior using the measurement from Aubourg et al. [42], Hq = 67.3 ± 1.1 km 
s -1 Mpc -1 . 


4 Results and discussion 


Since the SN la, BAO/CMB and H(z) data are independent from each other we can write 
the combined y 2 statistics as (see the Appendix) 


2 _ 2 ,2 . 2 
X — XSN ,marg ' XbAO/CMB XH.marg 


(4.1) 


where XsN,marg > X 2 bao/CMB and XH,marg are g ivei b respectively, by Eqs. (A.ll), (3.12) and 
(A. 20 ). At this point it should be pointed out that on any occasion when when we refer to 
an uniform prior marginalization over a parameter, we mean a marginalization over all its 
possible values. 

Let us first consider the case in which, besides fixing qi = 1/2, the final value of the 
deceleration parameter is fixed to the value qf = — 1 , such that asymptotically in the fu¬ 
ture (z = —1) a de Sitter phase is reached. Popular cosmological models like flat ACDM, 
quartessence Chaplygin, DGP etc, belong to this class of models as discussed in Ref. [4]. In 
Fig. 1 we show the 68 %, 95% and 99% confidence level regions in the (zt,r) plane. In all 
panels of this figure the blue dashed lines are contours for SN la + BAO/CMB (Planck) data. 
In the upper panel we also show the contours for SN la + BAO/CMB (Planck) + H(z) data, 
marginalizing over Hq with an uniform prior (purple lines). We see that including the H(z) 
data, assuming an uniform prior over iTo, does not change very much the contours nor the 
best fit. In the bottom panels we also show the results for the SN la + BAO/CMB (Planck) 
+ H(z) data, but marginalized over Hq with Gaussian priors. In the bottom left panel we 
used as prior for iTo, the Riess et al. [41] measurement (red lines) and in the bottom right 
we used the Aubourg et al. [42] result (green lines). Now the constraints are tighter, making 
clear the important role of the Hq prior when including H(z) data. We obtain that using 
the Riess et al result pushes the contours for lower values of r and higher values of zt when 
compared to the other two cases. We can quantify the improvement with respect to the pre¬ 
vious work [4] by comparing the areas inside the 95% C.L. contours. The areas ratio between 
the above results and the one from Giostri et al. [4], that used WMAP7, are 0.23, 0.23 and 
0.29 for the Riess, Aubourg and uniform priors. If we do not consider the H(z) data but use 
Planck instead of WMAP7 we get 0.31 for the areas ratio. It is also worth to mention that 
flat ACDM models (r = 1/3), represented by the black horizontal line in all figures, are in 
excellent accordance with the data in all cases analyzed. 

It is important to highlight the significant role of the new BAO data in the improvement 
of the constraints when considering the BAO/CMB test alone. The difference in the BAO 
dataset used in this work from that used in [4] (see Table 1 in this reference) is that here 
we dropped the measurements from [43] at z = 0.2 (which has one of the largest errors) and 
at z = 0.35. We considered instead the more precise estimates at z = 0.35 from [13] and at 
z = 0.57 from [14]. In Fig. 2 we compare our results for this test to the Giostri et al. [4] 
result, using CMB information from WMAP7 (the same used in [4]), WMAP9 and Planck 
2015 combined with the new BAO data. In the top left panel the only difference between 
the Giostri et al. result (blue dashed lines) and ours (purple lines) is the BAO data. We can 
see in the top right and bottom panels the same analysis replacing the WMAP7 information 
by WMAP9 (green lines) and Planck 2015 (red lines). We can quantify the improvement 
calculating the areas inside the contours. The 95% C.L. areas ratio relative to the Giostri 



*t z t 


Figure 1 . 68 %, 95% and 99% confidence contours in the plane z t vs. r, for qi — 0.5 and qf = —1. In 
all panels the blue dashed lines are the contours for SN la + BAO/CMB(Planck). In the top panel the 
purple lines are the contours for SN la + BAO/CMB (Planck) + H(z), marginalized over Hq with an 
uniform prior. In the bottom left panel the red lines are the contours for the same data, marginalized 
using the Riess et al. measurement as a Gaussian prior. In the bottom right panel the green lines are 
for the same data using the Aubourg et al. measurement as prior. In all panels the horizontal line 
corresponds to ACDM (r = 1/3) and the small dots and squares indicate the best-fit values. 


et al. result are 0.29 for Planck, 0.40 for WMAP9 and 0.47 for WMAP7. The numbers are 
similar for 68 % and 99% C.L.. Therefore, we conclude that the new BAO data gives a major 
contribution to the improvement found. This indicates that more precise BAO measurements 
are very important to achieve tighter constraints on the parameters. 

We show in Fig. 3 a comparison between the constraints for the case SN la + BAO/CMB 
+ H(z), marginalized over Hq with an uniform prior, using WMPA9 (green dashed lines) 
and Planck (blue lines). We can see that Planck’s higher precision gives us more stringent 
constraints. The 95% C.L. contour area ratio with Planck with respect to the one with 
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Figure 2. 68 %, 95% and 99% confidence contours in the plane Zt vs. r, for qi = 0.5 and qf = 
— 1, considering only BAO/CMB data. In all panels the blue dashed lines are the contours for 
BAO/CMB(WMAP7) data from Giostri et al. [4]. We show the results for our analysis, which uses 
a different BAO dataset, combined with CMB information from WMAP7 (in the top left panel, the 
purple lines), WMAP9 (in the top right panel the green lines) and Planck (in the bottom panel, the 
red lines). In all panels the horizontal line corresponds to ACDM (r = 1/3) 


WMAP9 is 0.71. We verified that this property also holds when the two other Hq priors are 
used. In the following we base our analyses by considering only the Planck information for 
CMB. 

We now consider the situation in which qf is free and can assume any value (we still fix 
qi = 1/2). In this very general case, the confidence regions of three model parameters (^, 
r and qf) have to be determined. However, following Ref. [4], instead of working with qf , 
we use go? the present value of the deceleration parameter. As discussed in Sec. 2 , Eq. (2.5) 
gives qf in terms of r and go allowing us to re-express q(z) in terms of these parameters 
(see Eq. (2.6)). It is important to mention that in the qf free case, to obtain the confidence 
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Figure 3. 68 %, 95% and 99% confidence contours in the plane z t vs. r, for qi = 0.5 and qf = — 1, 
using the combined data SN la + BAO/CMB + H(z), marginalized over Hq with an uniform prior. 
For the green dashed lines we used WMAP9 information while for the blue lines we used Planck 2015. 
The horizontal line corresponds to ACDM (r = 1/3) and the small dot and square indicate the best-fit 
values. 



uniform prior 

Aubourg et al. prior 

Riess et al. prior 


Qf = -1 

qf free 

Qf = - 1 

qf free 

Qf = - 1 

qf free 

Zt 

u.uu o Q3(o.07) 

n fi7 +0.10(0.22) 
u-u -0.08(0.13) 

n fi7 +0.03(0.05) 
u-u -0.03(0.06) 

n fi7 +0.10(0.22) 
-0.08(0.12) 

n fi o+0.03(0.05) 
u - u -0.03(0.06) 

PI 7f |+0.09(0.18) 
u ‘ ( u —0.08(0.13) 

T 

PI oq+0-04(0.09) 

^ -0.04(0.07) 

n 9(3+0-14(0.26) 
U ' ZU -0.10(0.18) 

PI q q - 1 - 0.04(0.08) 
0.03(0.06) 

n 9fi +0.14(0.26) 
U ' ZU -0.10(0.17) 

p, q.+0.03(0.07) 
U -' 3i -0.03(0.06) 

n „ .+0.13(0.25) 
U ‘ Z -0.09(0.16) 

qo 

- 

n .£+0.11(0.20) 
0.13(0.27) 

- 

n 4S +0.11(0.18) 
u '^°-0.13(0.26) 

- 

p| ,q + 0.11(0.19) 
L ’- OO -0.14(0.28) 


Table 1 . Summary of the best-fit values for all parameters when using SN la + BAO/CMB (Planck) 
+ H(z) , including the 68 % and 95% confidence intervals. For each parameter we marginalized over 
all other ones with an uniform prior assuming their all possible values. 


regions, following [4], we explicitly used the restriction on r given by Eq. (2.7). As discussed 
in Section 2, this restriction on r expresses the fact that the transition in the universe is from 
a decelerated phase to an accelerated one. 
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Figure 4. 68 %, 95% and 99% confidence contours in the plane r vs. z t for SN la + 

BAO/CMB(Planck) + H(z). The blue dashed lines in both panels are the contours for the case 
of marginalization over Hq with an uniform prior while the red ones (left panel) are for the Riess et 
al. prior and the green ones are for the Auborg et al prior. In all cases, the marginalization over 
the third model parameter is done with an uniform prior. In both panels the horizontal dashed line 
corresponds to ACDM (r = 1/3). 




Figure 5. 68 %, 95% and 99% confidence contours in the plane go vs. zt for SN la + 

BAO/CMB (Planck) + H(z). The lines and color conventions are the same described in Fig. 4 


In Figs. 4, 5 and 6 we display 2D marginalized (an uniform prior in the third free model 
parameter) confidence contours (68%, 95% and 99%) in the ( 2 ^, r), ( 2 ^, qo) and (go, t) planes 
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Figure 6. 68 %, 95% and 99% confidence contours in the plane go vs. r for SN la + 

BAO/CMB(Planck) + H(z). The lines and color conventions are the same described in Fig. 4 


for SN la + BAO/CMB (Planck) + H(z) data. In all panels of the three figures, the blue 
dashed lines are the contours for the case in which marginalization over Hq is performed with 
an uniform prior, while the red lines are the contours for the Riess et al prior and the green 
ones are for the Auborg et al prior. The 95% C.L areas ratio relative to the Giostri et al 
results are 0.37, 0.41 and 0.44, in the ( 2 ^, r) plane, 0.58, 0.58 and 0.71 in the ( 2 ^, go) plane 
and 0.67, 0.65 and 0.75, in the (go, r) plane, for the Riess, Aubourg and uniform priors, 
respectively. 

In Table 1 we show the best fit and the 68 % and 95% confidence limits on each free 
parameter after marginalizing over all the other parameters with an uniform prior. In the 
cases we fix qj — — 1 , there is no information in the Table on go since, in these cases, only 
zt and r are free parameters. However, we can use Eq. (2.5) to obtain the expected values 
of go- We obtain, with 68 % of confidence level, go = — 0 . 54 + 8 q 7 > #0 = —0.56^0 05 an d 
g 0 = -o.59l8;8« for the flat, Aubourg et al and Riess et al priors, respectively. 

In the left panel of Fig. 7 we show the r = 1/3 section of the 68 %, 95% and 99% confi¬ 
dence surfaces of the likelihood. We assume g/ is free and an uniform prior marginalization 
of H 0 . In the right panel, the r — 1/2 section is displayed. These two sections allow us 
to check directly the agreement between our results, for three free parameters ( 2 ^, r and 
go), and models with r fixed, like ACDM (r = 1/3 and gj = —1) and DGP (r = 1/2 and 
qj — —1) models, since the marginalized plots in Figs. 4, 5 and 6 can be misleading. We 
can see that ACDM, the black continuous line in the left panel, is in a very good agreement 
with the data, but the DGP model, the black dashed line in the right panel, is excluded at 
more than 99% confidence level. It is worth mentioning that the cut in the contours in the 
right panel is due to the constraint on the parameter r, given by Eq. (2.7), which ensures 
the transition from deceleration to acceleration, as discussed before. 

In Fig. 8 we show the r = 1/(3 x 0.785) (left panel) and r — 1/(3 x 1.216) (right panel) 
sections of the 68 %, 95% and 99% confidence surfaces of the likelihood. As in Fig. 7 we assume 
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Figure 7. Sections, at r = 1/3 (left panel) and r = 1/2 (right panel) of the 68%, 95% and 99% 
confidence surfaces in the parameter space (z t , r, go), marginalized over Ho with an uniform prior. 
The black continuous line (left panel) corresponds to ACDM and the black dashed line (right panel) 
corresponds to the DGP model. 




Figure 8. Sections r = 1/(3 x 0.785) (left panel) and r = 1/(3 x 1.216) (right panel) of the 68%, 
95% and 99% confidence surfaces in the parameter space (z t , r, go), marginalized over Ho with an 
uniform prior. The black continuous line in the left panel corresponds to reCDM with w = —0.785 
while that one in the the right panel corresponds to w = —1.216. 


that qf is free and an uniform prior marginalization over Hq. The black line in each panel 
corresponds to the flat w CDM model (w = p/p = constant), with w = —0.785 (left panel) and 
w = —1.216 (right panel). Flat teCDM are described by our parametrization if we assume 
r = — 4~. We can see from Fig. 8 that for wCDM our results indicate —1.22 < w < —0.78 at 
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more than 99% confidence level. 


5 Conclusions 

In this work we adopted a kinematical approach using a kink-like parametrization for the 
deceleration parameter to investigate the transition from cosmic deceleration to acceleration. 
To this end we considered recent SN la, BAO, CMB and H(z) observations. As in [1] and [4], 
we first considered the simple case in which qf — — 1 such that the models in this class have 
a final de Sitter attractor. With the new data the obtained constraints on the parameters 
are tighter but in good agreement with the previous ones. 

As compared with previous analysis [4] the first source for improvements is the SNela 
data used in this work. In [4] it was considered 288 SNela while in the present work we 
consider 740. A naive point of view could expect a significant improvement in the constraints 
due to the higher statistics. However, in the present work correlations and systematics 
uncertainties were more properly considered. It is important to highlight the impact of the 
systematic errors in the supernova data [12]. It is well known that systematics tend to 
raise the confidence contours and it is in fact currently dominant for SNela. Therefore, the 
SNela treatment considered in this work is more robust. As we mentioned in Sec. 4, we 
can get a 95% C.L. contour areas (for qj — —1) ratio, relative to the previous work, equal 
to 0.29 when using only BAO/CMB(Planck) data while we obtain 0.23 when we combine 
with SN la and H(z), which is an indication that the current SN la systematic errors are 
compensating most of the expected statistical enhancement. Regarding the BAO/CMB test 
we showed that, although we used more precise CMB information from Planck 2015 and 
WMAP9, the main source of improvements was the different BAO dataset used in this work. 
Here we disconsidered the older data at z = 0.2 and z = 0.35 and considered two more recent 
and precise measurements at z = 0.35 and z = 0.57. Of course, we expect to get tighter 
constraints once we have more and more precise BAO data in the future. 

With the inclusion of H(z) data we were able to study the impact of the prior on Hq 
in the results. We observed that assuming a Gaussian prior with the mean given by the 
measurement from Riess et al. favors slightly lower values of r and higher values of zt, 
as compared to the Aubourg et al. and the uniform prior marginalization. Although the 
ACDM model is in very good agreement with the results in all cases, the best fit in the 
Riess et al case is in slightly worse agreement. A possible explanation for this difference is 
that the Riess et al. measurement is based in low redshift observations and seems to be less 
model dependent than the Aubourg et al. one. Although, in large part of our analysis we 
adopted the more conservative approach based on an uniform prior marginalization over Hq, 
we note that imposing Gaussian priors on Hq is important to tighten the constraints on the 
parameters, mainly in the case in which qj = — 1. 

In the case of three free parameters ( 2 ^, r and qf) we also have improvements in the 
results as compared to previous work [4]. For instance, as can be seen in the right panel of 
Fig. 7, we are now able to exclude the DGP model at more than 99% of confidence level 
even considering systematic errors in the SN analysis. As illustrated in the left panel of 
Fig. 7, we observe once more that flat ACDM model is in very good agreement with all 
the observables considered in this work. Finally, we also showed that for w;CDM our results 
indicate —1.22 < w < —0.78 at more than 99% confidence level. 
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A Marginalization over the nuisance parameters 

In this appendix we describe how the marginalizations, mentioned before, are performed. 
We start with the full y 2 


x 2 — Xbao/cmb + Xsn + Xh> 


which corresponds to the likelihood 

T ( 1 2 

L oc exp I --x 


= exp \ -^Xbao/cmb ) ex P (~\xIn ) exp 


1 


\XOHD 


(A.l) 


(A.2) 

(A.3) 


that can be written as 


L(H 0 , M, 9) = L bao/CM b(0) L S n(Ho, M, 9) L h (H 0 , 9), (A.4) 

where 9 is a vector of the model parameters. In our case 6 — {z*,t, qf}. 

We now analytically marginalize over Hq and M such that 

L marg (9) = J L(H 0 ,M,9)dH 0 dM, (A.5) 

where the integral is considered over all the possible values of the parameters. The first 
difficulty we find, when integrating the expression above, is that L$n and Lh have a distinct 
functional dependence on Hq. The Xsn proportional to (logi^o) 2 while Xh proportional 
to Hq. Fortunately this problem can be solved by making the following change of variable: 
{Hq, M} —> {Hq, M = M + 5 log Ho}. So, the integral becomes 


Therefore, 


L 


marg 


( 0 ) 


L(H 0 ,M,9) 


d(H 0 ,M) 


dHo dM , 


JL SN (M,9 ) L h (H 0 ,9 ) L BAO/C MB{0)dH 0 dM , 

Lbao/cmb(6) J LsN(-Xi,9)dM J Ln(Ho,9)dHo, 

LBAO/CMB (^) LSN,marg{^B.marg (9). 


X 


2 

marg 


2 log Ljnarg ~ X^BAO/CMB + XSN ,marg + XOHD ,marg 


(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.10) 
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We now perform each marginalization as follows 


Lsn ,marg (0) ^ J ex P ipdsAj dM. oc / exp 


-- ( A-2MB + M 2 C) 


where 


A = (fji b - 51og D L (z b ,0)yC b - 51og D L (z b ,6)), 

li 


B = (fji b - 5 log D L (z b ,0)yC h 1, 

C = l t C 6 - 1 l. 

Solving the integral we get 

2 //-' 

XSN : marg ^ q ’ 

Similarly we obtain for Ln^marg 

L/H : marg OC J exp ^ dHo 


OC 


where 


/ exp -i (Aff 0 2 - 2i? 0 ^) 


^ _ ^2 E 2 ( z i, 0) 

iA a H,i 


dH 0 , 


N 


« = E 

2=1 


0)H o b s j 


a 


HA 


and E(z,6) — H t h(z,6 )/Hq. Thus, 


Xh 


,marg 


1 , B 2 
= log A — —. 


dM (A.11) 

(A.12) 

(A.13) 

(A.14) 

(A.15) 

(A.16) 
(A.17) 

(A.18) 

(A.19) 

(A.20) 


In the case we consider a Gaussian prior over Hq, we have to work with the posterior 


where 


P(H o ,M,0) = U{H o )L(H o ,M,0 ) , 

. Then the marginalized posterior is 

(A.21) 

J U(H 0 )L(H 0 ,M, 9)dH 0 dM 

(A.22) 

Lbao/cmb (0 ) J L sn {M, 9)dM J U(H Q )L H (H o ,0)dH Q 

(A.23) 

Lb AO/CM b(@) LsN,marg (^) PH.marg (^) • 

(A.24) 

PH,narg(0) = J P H (H 0 ,e)dH 0 

(A.25) 

= [ U(H o )L H (Ho,0)dHo 

(A.26) 
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and 


-2logP H (H 0 ,d) 


Xh 

N 


E 

2=1 

N 

E 

2=1 

N 


E 


-2iogn(ff 0 ) 

Ufh(Xi ! ffo, Vobs,i 

a H,i 

Hthjzu H [h 6) - Hpbsi 2 

Hthjzjr, H 0 ,9)- Hobsj, 2 


+ 


+ 




0 , 06 s 




Hthjz = 0,H 0 ,6) - Hobs ,0 
Tff,0 


The expression above has the same form as the xjj i n (3-18), but with one more data point, 
which corresponds to the prior. Therefore, we can proceed as before including only one more 
data point at z = 0. 
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